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Abstract 

We present a new analysis of the ratio e' /e which measures the 
direct CP violation m. K ^ vrvr decays. We use the l/Nc expan- 
sion within the framework of the effective chiral lagrangian for pseu- 
doscalar mesons. Prom general counting arguments we show that 
the matrix element of the operator is not protected from possi- 
ble large 'i/Nc corrections beyond its large- A'^c or VSA value. Cal- 
culating the l/Nc corrections, we explicitly find that they are large 
and positive. Our result indicates that a AI =1/2 enhancement is 
operative for Qg similar to the one of Qi and Q2 which dominate 
the CP conserving amplitude. This enhances e'/e and can bring the 
standard model prediction close to the measured value for central 
values of the parameters. 

1. Introduction. Recently, direct CP violation m K ^ tttt decays was 
observed by the KTeV |l[] and NA48 |0] collaborations. The new measure- 
ments are in agreement with the results of the NA31 experiment [Q. The 
present world average for the parameter e'/e is [0 

Re e'/e = (21.2 ±4.6) ■ 10-^ (1) 

which differs from zero by 4.6 standard deviations. In the standard model 
CP violation originates in the CKM phase, and direct CP violation is gov- 
erned by loop diagrams of the penguin type, in which the three quark 
generations are present. The value of e'/e is determined by a non-trivial 
interplay of the strong, electromagnetic, and weak interactions and depends 
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on almost all SM parameters. This makes its calculation quite complex. 
The main source of uncertainty in the calculation of the CP ratio is the 
QCD non-perturbative contribution related to the hadronic nature of the 
K 1T7C decays. Using the AS* = 1 effective hamiltonian, 

ntf'r = % ^« E Q^i^^) (/i < m,) , (2) 

the non-perturbative contribution, contained in the hadronic matrix el- 
ements of the four-quark operators Qi, can be separated, at the renor- 
malization scale ~ 1 GeV, from the perturbative Wilson coefficients 
Ci(/i) = Zi{fi) + Tyi{fi) (with r = -Xt/Xu and A, = V*^Vg^). Introduc- 
ing the matrix elements {Qi)i = {{mi)f\Qi\K), e'/e can be written as 

e' Gf oj IvdXt 

e 2 \e\ ReAn ' " ' ' ^ ' ' ^ ^ 



(3) 



u =Rey4o/Rey42 = 22.2 is the ratio of the CP conserving K ^ tttt isospin 
amplitudes; ^Iri+ri' takes into account the effect of the isospin breaking in the 
quark masses. In an exact realization of non-perturbative QCD, the scale 
dependences of the i/i and {Qi)i in Eq. must cancel, e'/e is dominated 
by (<56)o and {Q8)2 which cannot be fixed from the CP conserving data, 
different from most of the matrix elements of the current-current operators 
[|, ^. Beside the theoretical uncertainties coming from the calculation of 
the {Qi)i and of ^Ir^+r]', the analysis of the CP ratio suffers from the uncer- 
tainties on the values of various input parameters, in particular of the CKM 
phase in ImXt, of Aqcq = A^, and of the strange quark mass. The uncer- 
tainty on rris is especially important due to the density- density structure 
of Qq and Qs, which implies that their matrix elements are proportional to 
the square of the quark condensate and are hence proportional to I /ml. In 
all methods the calculation of the hadronic matrix element of Qq, which is 
a peculiar operator, appears to be the most difficult one. 

In this talk we present the results we obtained for e'/e together with 
CO. Kohler and E.A. Paschos We briefly explain the method used to 
compute the hadronic matrix elements, and discuss in detail the peculiari- 
ties in the calculation of the matrix element of Qq. 

2. General Framework. To calculate the hadronic matrix elements we 
start from the effective chiral lagrangian for pseudoscalar mesons which 
involves an expansion in momenta where terms up to 0{p'^) are included 0. 
Keeping only terms of 0{p'^) which contribute, at the order we calculate, 
to the K ^ TTTT amplitudes, for the lagrangian we obtain: 

Ceff = ^{{D,U^D^U) + ^^{\nU^-\nUf + {xU^ + Ux^)) 

+L,{D^U^D^U{x^U + U^x)) + L^{x^Ux^U + xU^U^) , (4) 
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with {A) denoting the trace of A, a = m^ + rn^i — 2m|., x = rA4, and = 
diag(m„, m^, rris) ■ f and r are parameters related to the pion decay constant 

and to the quark condensate, with r = — 2(gg)//^. The complex matrix 
?7 is a non-linear representation of the pseudoscalar meson nonet. The 
conventions and definitions we use are the same as those in @, ^, |^ • 

The method we use is the 1/Nc expansion introduced in |jlO|. In this ap- 
proach, we expand the matrix elements in powers of the external momenta 
and of l/Nc- From the lagrangian the mesonic representations of the quark 
currents and densities can be obtained by usual bosonization techniques. 
It is then straightforward to calculate the tree level (leading-iVc) matrix 
elements. For the l/N^ corrections to the matrix elements {Qi)i we cal- 
culated chiral loops as described in P, ^. Especially important to this 
analysis are the non-factorizable 1/Nc corrections, which are UV divergent 
and must be matched to the short-distance part. They are regularized by 
a finite cutoff which is identified with the short-distance renormalization 
scale [8-11]. The definition of the momenta in the loop diagrams, which 
are not momentum translation invariant, is discussed in detail in |^. 

For the Wilson coefficients we use both the leading logarithmic (LO) and 
the next-to- leading logarithmic (NLO) values Q]. In the pseudoscalar ap- 
proximation, the matching has to be done below 1 GeV. Values of the i/i at 
scales 0.6 GeV < /i < 0.9 GeV where communicated to us by M. Jamin. The 
NLO values are scheme dependent and are calculated within naive dimen- 
sional regularization (NDR) and in the 't Hooft-Veltman scheme (HV). The 
absence of any reference to the renormalization scheme in the low-energy 
calculation, at this stage, prevents a complete matching at the next-to- 



leading order ||12[. Nevertheless, a comparison of the numerical results 
obtained from the LO and NLO coefficients is useful as regards estimating 
the uncertainties and testing the validity of perturbation theory. 

3. Analysis of e'/e. At next-to-leading order, in the twofold expansion 
in powers of external momenta and of l/A^, we investigate the tree level 
contributions from the O(p^) and the C(p^) lagrangian as well as the one- 
loop contribution from the 0{p'^) lagrangian. For the matrix elements of 
the current-current operators, the corresponding terms are 0{p'^), O(p^), 
and 0{p'^/Nc), respectively; for density-density operators they are 0{p^), 
and C(p7AJ. 

Analytical formulas for all matrix elements at these orders are given 
in P, H. Among them four are particularly interesting and important: 



{Qi)o = ~^F^ {^K - 



4L^ o 1 



1 + -r^ml + 



- (An^F^ 



X 



6Ac-(^< + 6m2) logA^ + 



(5) 



3 



{Q 



2/0 



1 + 



4^5 2 



(4vr)2F2 

2 I 



(Q6)0 
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(6) 
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2V2 



-ml + logA^ + 



(4vr)2F2 



with R = 2m\ / {nis + rrid) . The elhpses denote finite terms which are 
not written exphcitly here. The constants L5 and Lg are renormahzed 
couphngs whose values are Lg = 2.07 ■ 10~^ and Lg = 1.09 ■ 10~^. Eqs. (|^) - 
(H) have several interesting properties H, First, the VSA values for 
{Qi)o and ((52)0 are far too small to account for the large AI = 1/2 en- 
hancement observed in the CP conserving amplitudes. Using the large- iVc 



limit [b[^^^^ = 3.05, B^^'''> = 1.22; 5)""'^ = Re{Qi)j/{Qi)J^^] improves 
the agreement between theory and experiment, but it still provides a gross 
underestimate. However, the non-factorizable 1/Nc corrections in Eqs. (P) 
and d^) contain quadratically divergent terms which are not suppressed 
with respect to the tree level contribution, since they bring about a factor 
of A = A^/ (47rF^)^ and have large prefactors. Varying Ac between 600 and 
900 MeV, b[^^^^ and B^2^^^ take the values 8.2 - 14.2 and 2.9-4.6, respec- 
tively. Quadratic terms in {Qi)o and {Q2)o produce a large enhancement 
which brings the AI =1/2 amplitude in agreement with the data ||^. Cor- 
rections beyond the chiral limit {rriq = 0) in Eqs. (|^) - (||) are suppressed by 
a factor of 6 = rnj^^/ (47rF^)^ ^ 20 % and were found to be small. 

For the operators Qq and Qs, values rather close to the VSA [Bq'^'^^ = 
Bf'^ = 1] are used in the literature. As a result the experimental range for 
e' / e can be accommodated in the standard model only if there is a conspir- 
acy of the input parameters m^, flr^+r]', ImAj, and Aqcd (see e.g. f^). The 
fact that the VSA fails completely in explaining the AI = 1/2 rule therefore 
raises the question whether it can be used for Qq and Qs- In fact, at the 
present stage of the calculation, the case of {Qq)o and {Q8)2 is different from 
that of (Qi,2)o- The leading- Ai'c values are very close to the corresponding 
VSA values. Moreover, the non-factorizable loop corrections in Eqs. (|^) 
and d^), which are of 0{p^/Nc), are found to be only logarithmically diver- 
gent p|. Consequently, in the case of {Q8)2 they are suppressed by a factor 
of 6 compared to the leading 0{p^) term and are expected to be of the order 
of 20 % to 50 % depending on the prefactors. We note that Eq. (H) is a full 



>{l/2) 



?(A/) 
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leading plus next-to- leading order analysis of the Qg matrix element. The 
case of B^^'^^ is more complicated since the 0{p^) term vanishes for Qg- 
Nevertheless, the non-fact or izable loop corrections to this term remain and 
have to be matched to the short- distance part of the amplitudes These 
0{p^/Nc) non-factorizable corrections must be considered at the same level, 
in the twofold expansion, as the 0{p'^) tree contribution. Consequently, a 
value of B^^^"^^ around one [which corresponds to the 0{p^) term alone] is 
not a priori expected. However, numerically it turns out that the 0{p^/Nc) 
contribution is only moderate. This property can be understood from the 
{U'')dq{U)gs structure of the Qq operator which vanishes to 0{p^) implying 
that the factorizable and non-factorizable 0{p^ /Nc) contributions cancel to 
a large extent P]. This explains why the deviations we observe with respect 
to the VSA values are smaller than for Qi and Q2 and why in particular 
for Qq to 0{p^/Nc) we do not observe a AI =1/2 enhancement. Varying 

Ac between 600 and 900 MeV, B^^^^^ and 5^ take the values 1.10 - 0.72 

(1/2) (3/2) 

and 0.64 — 0.42, respectively. Bq and Bq are therefore more efficiently 
protected from possible large 1/Nc corrections of the O(p^) lagrangian than 
bI^C^K The effect of the 0{p^/Nc) term is however important for i?g^^^^ as 

(3/2) 

for Bq because it gives rise to a noticeable dependence on the cutoff 
scale . We note that B^^^"^^ shows a scale dependence which is very simi- 
lar to the one of Bq^^"^^ leading to a stable ratio Bq'^'^^ / B^^"^^ ~ 1.72 for Ac 
between 600 and 900 MeV. The 0{p^/N^ corrections consequently make 
the cancellation of and Qs in ^'/^ l^ss effective, but the values of the 
matrix elements are reduced. Hence we obtain values for e' /e @ close to 

/I/O) ('3/2) 

the ones found with the VSA values of and -Bg ; in particular, for 

central values of the input parameters [ms(l GeV)=150 MeV, Vt^+j^i = 0.25, 
ImAt = 1.33 ■ 10""^, and Aqcd = 325 MeV, see P and references therein] the 
values for the CP ratio are significantly smaller than the data. 

However, since the leading 0{p^) contribution vanishes for Qq, correc- 
tions from higher order terms beyond the 0{p^) and 0{p^/Nc) are expected 
to be large 0. The terms of O(p^) and 0{p^/Nc) correspond to the lowest 
(non- vanishing) order, and the calculation of the next order terms is very 
desirable. In the twofold expansion, the higher order corrections to the ma- 
trix element of Qe are of orders: 0{p^), 0{p°/N^), and 0{p^/Nc). A full 
calculation of these terms is beyond the scope of our study. In particular, 
higher order terms in the p^ expansion, which are chirally suppressed, can- 
not be calculated because the low-energy couplings in the 0{p^) lagrangian 
are very uncertain or even unknown. In we investigated the 0{p'^/Nc) 
contribution, i.e., the l/N^ correction at the next order in the chiral expan- 
sion, because it brings about, for the first time, quadratic corrections on the 
cutoff. We remind the readers that for the CP conserving amplitude it is 
mainly the (quadratic) OijP' /N^ corrections which bring to (Qi,2)o a large 
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enhancement relative to the (leading-A^^c) values. As the leading-A'^c 

value for Qq is also O(p^) we cannot a priori exclude that the value of {Qq)o 
is largely affected by 0{p'^/Nc) corrections too, removing from Qq the prop- 
erty observed to 0{p^/Nc), to be protected from large 1/Nc corrections. As 
explained above, quadratic 0{p'^/Nc) corrections are proportional to the 
factor A = A^/(47rF^)^ relative to the 0{p'^) tree level contribution. Dif- 
ferent is the case of since its leading-iVc value is 0{p^) at lowest order 
in the chiral expansion. Quadratic terms for Qs are consequently chirally 
suppressed with respect to the leading- A^c value. 

Calculating the quadratic term of 0{p'^/Nc) for matrix element of Qq 
and adding it to the and OijP /Nc) result in Eq. (|^) we obtain: 



(Q6)C 



4v^ 



K 



\ 3 logAg 

2 (477)2^2 ) 16 (477)2 



+ 



(9) 



The result for the 0{p'^/Nc) term we already presented in |T^. Numerically, 
we observe a large positive correction from the quadratic term in Eq. @. 
The slope of this correction is qualitatively consistent and welcome since it 
compensates for the logarithmic decrease at (9(p°/Ac). Varying A^ between 
600 and 900 MeV, the B^^^"^^ factor takes the values 1.50 — 1.62. The approx- 
imate stability of B^^^"^^ is in accordance with the perturbative evolution, 
since the non-diagonal dependence of uq on the renormalization scale, i.e., 
the one beyond the leading- A'c scale dependence of E? in Eq. (^, was found 
to be small. The quadratic term of 0{p'^/Nc) is of the same magnitude as 
the 0{p'^) tree term. Qg is a A/ = 1/2 operator, and the enhancement of 
{Q%)q indicates that at the level of the 1/Nc corrections the dynamics of 
the A/ = 1/2 rule applies to as to Qi^2- One might however note that 
the enhancement observed for Qq is smaller than for Qi and Q2- 

Using the quoted values for Bq^^'^^ together with the full leading plus 
next-to-leading order B factors for the remaining operators we calculated 
e' /e for central values of m^, firj+r?', IniAt, and Aqcd- The results for the 
three sets of Wilson coefficients LO, NDR, and HV and for Ac between 
600 and 900 MeV are given in Tab. |I]. The numbers are obtained with two 
different methods for analyzing the sensitivity on the imaginary part coming 
from the final states interactions. In the first case, we use the real part of 
our calculation and the phenomenological phases 5o = (34.2±2.2)° and 82 = 
(-6.9±0.2)° |ll,andreplace|E^2/^(Q^)/|inEq. (|) by E^ l/iRe(Qi)// cos 5/ 
[ |15|1 . In the second case, we use only the real part assuming zero phases. 
The latter case is very close to the results we would get if we used the small 
imaginary part obtained at the one- loop level 0. Collecting together the 
LO, NDR, and HV results for the two cases and for central values of the 
parameters, we find the following (conservative) range: 



7.0 -lO"^ < e'/e (central) < 24.7-10" 
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Case 1 


Case 2 


LO 
NDR 
HV 


19.5 < e'je < 24.7 
16.1 < e'/e < 23.4 
9.3 < e'/e < 19.3 


14.8 < e'/e < 19.4 
12.5 < e'/e < 18.3 
7.0 < e'/e < 14.9 


LO + NDR + HV 


9.3 < e'/e < 24.7 


7.0 < e'/e < 19.4 



Table 1: Central ranges for e'/e (in units of 10 ^) as explained in the text. 



which is in the ball park of the experimental result in Eq. (|T]). Performing 
a complete scanning of the parameters [125 MeV < ms{l GeV) < 175 MeV, 
0.15 < 1]^+^/ < 0.35, 1.04 ■ 10"^ < ImAt < 1.63 ■ 10-^ and 245 MeV < 
Aqcd < 405 MeV] we obtain 2.2 • 10"^ < e'/e (scanned) < 63.2 ■ 10"^ (see 
Tab. 1^). The numerical values in the tables can be compared with the re- 
sults of the Munich, Trieste, and Rome groups 0, |15|, |16|]. The values for 
B^^^"^^ can also be compared with |]1^. The large ranges reported in Tab. |] 



can be traced back to the large ranges of the input parameters. The param- 
eters, to a large extent, act multiplicatively, and the larger range for e'/e is 
due to the fact that the central value (s) for the ratio are enhanced roughly 
by a factor of two compared to the results obtained with B factors for Qq 
and Qs close to the VSA. More accurate information on the parameters, 
from theory and experiment, will restrict the values for the CP ratio. 

4. Summary. We have shown that the operator Qq, similar to Qi and Q2, 
is not protected from large 1/Nc corrections coming from quadratic terms 
of 0{p'^/Nc). From general counting arguments we have good indications 
that among the various next-to- leading order terms in the and 1/Nc 
expansions they are the dominant ones. Calculating those terms we find 
that they enhance B^^'^^ and bring e' / e much closer to the data for central 
values of the parameters. We obtain a quadratic evolution for Qq which 
indicates that a A/ = 1/2 enhancement is operative for similar to the 
one of Qi and Q2- Bf^"^^ is expected to be affected much less by terms 
of 0{p^ /Nc) due to an extra suppression factor relative to the leading 
(9(p°) tree term. 

One should recall that our analysis of the 0{p'^/Nc) terms for Qq is 
performed in the chiral limit. It would be desirable to calculate the correc- 
tions beyond the chiral limit, from logarithms and finite terms. It would 
also be interesting to investigate the effect of higher resonances. Each of 
the additional effects separately is not expected to counteract largely the 
enhancement found for B^^'^\ Nevertheless, in the extreme (and unlikely) 
case where all these effects would come with the same sign a significant mod- 
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Case 1 


Case 2 


LO 
NDR 
HV 


8.0 < e'je < 62.1 
6.8 < e'je < 63.9 
2.8 < e'le < 49.8 


6.1 < e'/e < 48.5 

5.2 < e'/e < 49.8 
2.2 < e'/e < 38.5 


LO + NDR + HV 


2.8 < e'/e < 63.9 


2.2 < e'/e < 49.8 



Table 2: Same as in Tab. but for the complete scanning of the parameters. 



ification of the result cannot be excluded formally. In order to reduce the 
scheme dependence in the result for e'/e, appropriate subtractions would 
be necessary (see [|17|, |T8|). 
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